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Abstract 



We simulated the million degree interstellar medium 
and its soft X-ray images in the disk and halo of spiral 
galaxies using the bipolar hypershell (BHS) model. In 
this model dumbbell- or hourglass-shaped expanding 
shells of several kpc radii are produced by a sudden 
energy release in the central region. We then applied 
our model to a mini-sample of starburst galaxies seen 
under different inclinations, namely the nearly edge-on 
galaxies NGC 253, NGC 3079 and M 82, the highly 
inclined galaxy NGC 4258 as well as the nearly face- 
on galaxies NGC 1808 and M 83. For all galaxies, our 
results reproduce the X-ray characteristics observed in 
the 0.1-2.4 keV ROSAT energy band: The bipolar hy- 
pershell morphology, the spectral energy distribution 
of the diffuse disk and halo emission as well as absorp- 
tion gaps in the diffuse X-ray emission caused by a 
shadowing of soft X-rays due to cold intervening gas in 
the disks of the galaxies. In general, the required total 
energy for the starburst is estimated to be of the order 
of 10 55 ergs, corresponding to the overall kinetic en- 
ergy generated in ~ 10 4 type-II supernova explosions. 
The expansion velocity of the shells is estimated to be 
~ 200 km s -1 , which is necessary to heat the gas to 
~ 0.2 keV (2.3 million K), and the age to be of the 
order of 3 x 10 7 years. In the case of the very nearby, 
nearly edge-on galaxy NGC 253 all characteristics of 
the BHS model can be studied with high spatial reso- 
lution. 

Using the property that the shell morphology is sen- 
sitive to the ambient density distribution, we propose 
using soft X-ray data to probe the gas distributions in 
the disk, halo and intergalactic space in general. The 
application of our model to images at higher spatial and 
spectral resolution, as provided by Chandra and XMM, 
will help us to further disentangle the ISM density dis- 
tributions and will lead to a better understanding of 
the disk halo interface. 
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I Introduction 

A bipolar hypershell (BHS) model has been proposed 
to explain the double-horn features observed in the 
radio-continuum halo of the starburst galaxy NGC 253 
(Sofue 1984). Recent soft-X-ray observations of NGC 
253 with ROSAT revealed extended outer-halo emis- 
sion with double-horn morphology, and is interpreted 
as being a hypershell driven by a central starburst 
(Vogler and Pietsch 1999a; Pietsch et al. 2000). Many 
spiral galaxies exhibit similar galactic-scale outflows 
in the form of dumbbell-shaped hypershells, bipolar 
cylinders, or galactic-scale jets. Examples of proto- 
types of such galaxies are M 82 (NGC 3034), NGC 
1808 and NGC 3079 (see next section for references). 
Starburst-origin models for extended halo features in 
galaxies were proposed by several authors (Heckman 
et al. 1993; Tomisaka and Ikeuchi 1987; Suchkov et al. 
1994, 1996). 

The Milky Way also exhibits bipolar hypershells ex- 
tending above and below the galactic plane in soft X- 
rays (Snowden et al. 1997) and radio continuum (e.g. 
Haslam et al. 1982), which are also interpreted as an 
effect of a giant Jl-shaped shock front originating from 
a starburst some 10 7 years ago in the Galactic Center 
(Sofue 1977, 1984, 1994, 2000). These hypershells and 
flows, including that in the Milky Way, require a similar 
amount of total energy of the order of ~ 10 55 ergs cor- 
responding to the energy release from 10 4 to 10 5 type 

II supernovae. This fact suggests that they may have 
a common property and origin, and can be modeled by 



1 



a unified scheme of energy release in the galactic cen- 
ter and its propagation into the halo and intergalactic 
space. 

In this paper, we simulate the extended soft X-ray 
halo components in various galaxies based on the bipo- 
lar hypershcll (BHS) model developed for the Milky 
Way halo (Sofue 1977, 1984, 1994, 2000). We further 
propose to use the simulation to date and measure the 
explosive events at the nuclei, such as their starburst 
history. The simulation can also be used to obtain 
information about the galactic-halo-intergalactic inter- 
face in disk galaxies. The X-ray absorption morphol- 
ogy near the galactic plane can be used to probe the 
interstellar gas distribution in the galactic disk. 

2 Soft X-ray Hypershells in 
Nearby Spiral Galaxies 

2.1 NGC 253 

ROSAT images of the nearby (d — 2.6 kpc, Puche & 
Carignan 1988) starburst galaxy NGC 253 revealed 
diffuse X-ray emission from the nuclear region, the 
disk and the halo hemispheres (Pietsch et al. 2000). 
The emission measures of the three components are 
nearly identical, and the integral diffuse emission (L x = 
4 x 10 39 erg s -1 ) accounts for <~ 80% of the total X-ray 
luminosity from NGC 253. The diffuse bulge emission 
is due to an extended, highly absorbed source with an 
extent of 250 kpc, and a hollow-cone shaped "X-ray 
plume" with an extent of ~ 700 pc along the SE mi- 
nor axis. The X-rays in the bulge region are thought 
to trace the interaction between a galactic super-wind 
driven by the starburst nucleus and dense interstellar 
medium in the disk. 

In the disk, the highly absorbed diffuse emission fol- 
lows the spiral arm structure. It most probably re- 
flects a superposition of diffuse X-ray emission from 
hot interstellar and discrete sources, such as emission 
from binaries, SN-remnants or Hii regions. Two major 
components of the diffuse halo emission can be sepa- 
rated: "coronal emission" arising from hot interstel- 
lar medium close to the disk (scale height ~ 1 kpc) 
and emission from the "outer halo" . NGC 253 is not 
completely oriented edge-on, and the NW side of the 
galaxy is closer to us. Due to this viewing geome- 
try and absorption of the soft X-rays by the disk of 
NGC 253, the coronal emission can only be detected 
from the near (SE) side of the galaxy, while emission 
from the outer halo component is partly shielded by the 
Hi disk of NGC 253 to the NW. Following Pietsch et 
al. (2000), the X-ray corona might be fueled by hot gas 
from galactic fountains within the disk, while the outer 



halo emission is attributed to a superwind emanating 
from the starburst nucleus. The latter assumption can 
also explain the "hour-glass-" or "dumbbell-shaped" 
morphology of the outer X-ray halo. Pietsch & Voglcr 
(2000) estimate temperature of the outer halo compo- 
nent at 0.2 keV (2.3 million K), corresponding to an 
equivalent shock- heating velocity of 200 km s _1 . 

2.2 NGC 1808 

NGC 1808 is a galaxy with large-scale cylindrical dusty 
jets emanating from the nucleus at high velocities, and 
is supposed to be a mildly tilted M82 type galaxy 
(Veron-Cetty and Veron 1985; Phillips 1993). Soft X- 
ray imaging observations with the ROSAT PSPC and 
HRI revealed an active X-ray nucleus and extended 
emission in the disk region (Junkes et al. 1995; Dahlem 
et al. 1994), as well as a very extended X-ray emission 
in the halo (Iwaki et al. 1996). The total soft X-ray 
luminosity is of the order of 3 x 10 40 ergs s _1 . Iwaki 
et al. (1996) obtained wide area ASCA images in the 
0.5-2 keV band, and found that the extended compo- 
nent is highly elongated in the direction of the dusty 
jets on both sides of the galactic plane. Furthermore, 
the X-ray halo appears to be bifurcated into two ridges, 
suggesting that it comprises double-horn feature or a 
partial shell structure. 

2.3 NGC 3079 

The nearly edge-on galaxy NGC 3079 is known to host 
an active nucleus classified MS iX LINER or Seyfert n 
nucleus, and giant radio lobes and Ha superbubbles 
emanate along the minor axis (Irwin & Seaquist 1991; 
Cecil et al. 1995). ROSAT observations (Pietsch et al. 
1998) detected an integral 0.1-2.4 keV luminosity of 
L x = 3 x 10 40 erg s _1 , which can be resolved into three 
components. One component consists of extended 
emission from the central region (L x = 1 x 10 40 erg s _1 ) 
originating from the region of the central superbubble. 
While most of this emission is thought to represent hot 
gas in the superbubble, one cannot exclude additional 
contributions of < 30% from the active nucleus. The 
second main emission component arises from the disk 
of the galaxy (L x = 7 x 10 39 erg s _1 ) and can be partly 
resolved into point-like sources. Each of the three de- 
tected point sources has L x ~ 6 x 10 38 erg s _1 . Finally, 
emission from the halo (i x = 6 x 10 39 erg s _1 ) is de- 
tected. As in the case of NGC 253, the hot gas does not 
fill the halo hemispheres uniformly, but an hourglass- 
shaped structure of the halo emission is reported. A 
mean temperature for the halo gas was estimated to be 
3.5 x 10 6 K. Interestingly, the lowest contours of the 
extended emission show a saddle shape, with a neck at 
the galactic plane, suggestive of bipolar double-horns. 
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2.4 NGC 4258 

NGC 4258 's anomalous arms could also be due to 
out-of-planc partial shells (van Albada and van der 
Hulst 1982; Vogler and Pictsch 1999b). Diffuse X- 
ray emission of the galaxy was detected with L 
2 x 10 40 erg s" 1 and fills - 40% of the D 25 diameter 
of NGC 4258. At least one half of the diffuse emission 
is caused by the anomalous spiral arms of this galaxy, 
and the X-ray emission morphology closely follows the 
pattern detected in the radio wavelengths. 

In the literature, those arms are discussed to be (a) 
jets confined to the disk of the galaxy and expelled by 
the active nucleus (e.g., Martin ct al. 1989, Hummel 
et al. 1989, Plante et al. 1991, Cecil ct al. 1995), (b) 
outflows (jets or plumes) from the active nucleus leav- 
ing the plane of the galaxy and possibly 'raining down' 
onto the disk again (van der Kruit et al. 1972), (c) a 
hypershell of hot interstellar medium as present in the 
nuclear region of NGC 253 and expanding from the 
disk into the halo (Sofue 1984, 2000), or (d) bar shocks 
heating the interstellar medium due to the high differ- 
ential velocities between the bar orbit and interstellar 
medium in the bulge and inner disk of NGC 4258 (e.g., 
Cox & Downcs 1996). 

Remaining diffuse emission not correlated with the 
anomalous arms might originate from the disk and the 
halo of the galaxy. For a possible halo component, 
the temperature was estimated to be 2 x 10 6 K. This 
temperature is lower by a factor of 2 than the temper- 
ature of the hot gas along the anomalous spiral arms 
(T > 4 x 10 6 K). 

In addition to the diffuse emission, 14 point sources 
(integral L x ~ 2 x 10 39 erg s" 1 ) were found in the 
galaxy. For possible contributions of the active nu- 
cleus to the X-ray emission, an upper limit of L x ~ 
1 x 10 38 erg s- 1 was found in the ROSAT 0.1-2.4 keV 
band. 

2.5 M 82 

M82 is an edge-on starburst galaxy, which ejects 
a galactic-scale flow through bipolar cylindrical jets 
(Nakai et al. 1987) to which larger-scale X-ray ex- 
tended features are associated (Lehnert et al. 1999). 
Strickland et al. (1997) report the detection of a galac- 
tic wind. The X-ray emission was visible out > 6 kpc 
above and below the plane of the galaxy, and the lumi- 
nosity of the diffuse X-ray emission outside the nuclear 
region is ~ 2 x 10 40 erg s _1 . From X-ray spectra, the 
temperature of the gas is found to vary from 0.6 keV 
near the nucleus to 0.4 keV for distances further away 
from the plane. The authors compared the X-ray re- 
sults with several theoretical models, namely Tomisaka 
& Ikeuchi (1988), Tomisaka & Bregman (1993) and 



Suchkov et al. (1994). From their comparison, they 
favor the model of Suchkov et al. (1994). This model 
assumes a two component cold rotating dense disk, a 
non-rotating hot tenuous halo, and a starburst history 
incorporating the milder mass and energy input from 
stellar winds into the halo before more energetic super- 
nova explosions dominate the scenario. Shocked halo 
gas then provides the major emission component in the 
ROSAT band. 

A reanalysis of the ROSAT data (Lehnert et al. 
1999) and new Ha images showed the existence of a 
4 kpc x 1 kpc region of spatially coincident X-ray and 
Ha emission about 10 kpc north of the nucleus. The 
luminosity in both wavebands accounts for roughly 1% 
of the total M 82 luminosities in both bands. The au- 
thors attribute the Ha emission to ionizing radiation 
from the starburst that propagates into the halo, and 
the X-ray emission to a shock-heated massive ionized 
cloud in the halo, which is encountered by the super- 
wind. With this scenario, they follow theoretical pre- 
dictions of Suchkov et al. (1996). 



2.6 The Milky Way 

Using ROSAT All Sky Survey (RASS) data, presented 
in several energy bands from 0.1 to 2.4 keV in Snow- 
den et al. (1997), we showed that the North Polar 
Spur (NPS) and its western and southern counter-spurs 
trace a giant dumbbell-shape wit a neck at the galac- 
tic plane (Sofue 2000). These features are interpreted 
as effects of a shock front originating from a starburst 
15 million years ago with a total released energy of 
the order of ~ 10 56 ergs or 10 5 type II supcrnovae. 
We simulated all-sky distributions of radio continuum 
and soft X-ray intensities based on the bipolar hyper- 
shell galactic center starburst model. The simulations 
nicely reproduce the radio NPS and related spurs, as 
well as radio spurs in the tangential directions of spiral 
arms. Simulated X-ray maps in 0.25, 0.75 and 1.5 keV 
bands reproduce the ROSAT X-ray NPS, its western 
and southern counter-spurs, and the absorption layer 
along the galactic plane. We further proposed to use 
the ROSAT all-sky maps to probe the physics of gas in 
the halo-intergalactic interface of the Milky Way, and 
to directly date and measure the energy of a recent 
Galactic Center starburst. The present paper is a se- 
ries of the BHS modeling of extended halo X-rays in 
galaxies using the same numerical code. 
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3 Bipolar Hypershell Starburst 
Model 

We adopt the same numerical methods to trace the 
shock envelope and X-ray intensity distributions as we 
used to model the radio spurs and BHS in the Milky 
Way, which is described in detail in Sofue (2000). We 
briefly describe the methods below. 

3.1 Adiabatic-Shock Envelope Method 

The propagation of a shock wave through the galac- 
tic halo induced by a point energy injection at the 
center can be calculated by applying the shock enve- 
lope tracing method of Sakashita (1971) and Mollcnhoff 
(1976), who extended the Laumbach and Probstein 
(1969) method for tracing the evolution of a shock front 
to a case of axi-symmetric distributions of ambient gas. 
The flow field is assumed to be locally radial, and the 
gas is adiabatic, and, therefore, the heat transfer by 
radiation and counter-pressure are neglected. The den- 
sity contrast between the shock front and ambient gas 
is given by (7 + l)/(7 — 1) = 4 for 7 = 5/3, where 7 is 
the adiabatic exponent of the gas. 

The equation of motion of the shock wave is given 
as follows (Mollenhoff 1976): 

f R P f R 1 f dr\ 2 

E = l —^dr + J o Po47 n> (1). 

Here, E is the total energy of the explosive event, P 
is the internal pressure, 7 is the adiabatic exponent, 
po is the unperturbed ambient gas density, and r is the 
radius from the explosion center, with suffix denoting 
the quantities of the unperturbed ambient gas, and R 
is the radius of the shock front. This equation then 
leads to an equation of the shock radius, as described 
in Sofue (2000). 

The unperturbed density distribution of gas in a 
galaxy is assumed to comprise a stratified disk, a halo 
with an exponentially decreasing density, and inter- 
galactic gas with a uniform density. The differential 
distribution is approximated by the following expres- 
sion. 

pa = picxp(-(z/zi)' ) )) + p 2 cxp(-z/z 2 ) + p 3 (2). 

Here, suffices 1, 2 and 3 denote quantities for the disk, 
halo and intergalactic gas, respectively, p is the den- 
sity, z is the height from the galactic plane, Zi is the 
scale thickness of the disk and halo. The power index 77 
represents the tightness of the disk component toward 
the galactic plane, and is taken either 1 (exponential) 
or 2 (gaussian). The third component, representing 
the intergalactic gas density, is taken to be constant 



at ps = 10 _5 H cm -3 , while the other parameters are 
set in various possible combinations as listed in Table 
1. We also calculate a case in which the halo density 
decreases with a power law. We further calculate an- 
other case with uniform halo density which decreases 
suddenly at a certain height facing the intergalactic 
gas, which represents a case for halo gas confined by 
the intergalactic gas pressure. 

Solutions of the above equations are presented in Fig. 
1 for various parameter sets given in Table 1. The 
total energy input is taken to be E n = 0.75 x 10 55 
ergs for SI to S5, and 3 x 10 55 ergs for S6, and the 
shock envelopes are drawn every 10 million years. The 
shock velocity can be estimated by dividing the radial 
distance between succeeding envelopes by 10 7 years. 

Models SI and S2 : Exponential disk with scale 
height 0.1 and 0.2 kpc, respectively, exponential halo 
with scale height 1 kpc, and intergalactic gas of uni- 
form density. 

Models S3 and S4 : Gaussian disk with half thickness 
0.1 and 0.2 kpc, respectively, exponential disk of scale 
height 1 kpc, and intergalactic gas of uniform density. 

Model S5: Power-law density disk, power-law den- 
sity halo, and intergalactic gas of uniform density. 

Model S6: Disk of 2 kpc thickness with sudden de- 
crease of density, facing the intergalactic low-density 
uniform gas. This model represents a case where the 
halo-to-disk gas is not gravitationally stratified (e.g., 
as exponentially), but is confined by intergalactic gas 
pressure, so that the internal pressure and density in 
the halo are more uniform. 
- Fig. 1 - 

3.2 X-ray Emission and Absorption 

The X-ray intensity distribution for a model galaxy 
having a hypershell is calculated as below. The centers 
of the shells are assumed to be at z — ±6 kpc on the 
rotation axis, and the radii 6 and 9 kpc in the radial 
and vertical directions, respectively. The volume emis- 
sivity is calculated from the density contrast of the 
shocked gas, which is assumed to be proportional to 
p/po = (7+l)/(7 — 1), where p and po are gas densities 
in the shocked shell and unperturbed halo gas, respec- 
tively. In the present simulation, the profile of cmissiv- 
ity perpendicular to the shell surface is simply repre- 
sented by an exponentially decreasing function behind 
the shock front toward the center with a scale thick- 
ness of 500 pc. The emissivity also decreases with the 
height from the galactic plane with a scale height of 3 
kpc, corresponding to the exponentially decreasing gas 
density in the halo. 

The X-ray emission from a hypershell is assumed to 
be thermal free-free radiation, whose emissivity is given 
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Table 1: Parameters for Shock Envelope Solutions 



Model 


Energy 


Pi 


hi 


Disk z type 


92 


h 2 


P3 




ergs 


H cm" 3 


kpc 




H cm" 3 


kpc 


H cm- 3 


SI 


0.4 x 10 55 


1 


0.1 


exponential 


0.01 


1 


10" 5 


S2 


11 


1 


0.2 


exponential 


0.01 


1 


11 


S3 


11 


1 


0.1 


gaussian 


0.01 


1 


11 


S4 


11 


1 


0.2 


gaussian 


0.01 


1 


11 


S5 


11 


1 


0.1 


power law 


0.01 


1 


11 


S6 


1.6 x 10 55 


0.1 


2 


exph^/M 4 ] 


0.0 




11 



The shock envelope (radius and velocity) obeys a scaling law: For equal ratios E/pi, one obtains the same 
envelope, independent of the individual E, pi values 



by 

e = n 2 A oc p 2 T^ 2 (3). 

Here, A is the cooling function of the gas. Since the 
temperature is as high as T ~ 10 7 K in our simulation, 
the radiation is almost totally free-free, and the contri- 
bution by recombination lines from hydrogen, helium, 
and metals, which are significant at T ~ 10 5 K, can be 
neglected. A typical emission measure along the hy- 
pershell ridge is n 2 L ~ 10~ 2 cm~ 6 pc for a tangential 
pass in the shell of about 2 kpc. For simplicity, we 
assume that the temperature of X-ray emitting gas is 
constant in the shocked shell. In addition to the hy- 
pershells, we assume a galactic disk X-ray component 
of scale height of 500 pc, and a bulge component of 
scale radius 1 kpc. X-ray emissions from the ambi- 
ent (unperturbed) halo gas is neglected, because the 
density would not be sufficiently high to contribute a 
significant 0.1-2.4 keV X-ray emission measure. 
The X-ray intensity, /, is calculated by 

dl = exds — Kids (4), 

where the first term in the right-hand side is the emis- 
sion measure with ex being the X-ray cmissivity, and 
the second term represents the absorption rate with k 
being absorption coefficient, and s is the distance along 
the line-of-sight. 

The interstellar extinction of soft X-rays occurs due 
to the photoelectric absorption by metals, and has been 
calculated for the solar metal abundance by Morrison 
and McCammon (1983). The absorption coefficient k 
is, then, represented by 

K =(n H /A£)(£ X /£ X )- 2 - 5 (5), 

where nu is the number density of hydrogen atoms, and 
is the e-folding column density of interstellar neu- 
tral hydrogen at the photon energy i? x (Ryter 1996; 
Sofue 2000). At Ef = 0.75 keV, the e-folding column 
density for Ef = 0.75 keV is given as JVg = 3 x 10 21 



H cm~ 2 (Snowden et al. 1997). Since the ROSAT en- 
ergy bands are rather broad, covering higher and lower 
energies around the representative energies, we adopt 
two representative cross sections at 0.75 and 1.5 keV 
bands, which are 1 and 0.1 times the value at 0.75 keV, 
respectively. 

The absorbing disk gas is assumed to be condensed 
in logarithmic spiral arms of pitch angle 6° and the 
width is taken to be 1/5 of the arm separation, as in 
the same way assumed for the Milky Way. The peak 
hydrogen density in an arm at a distance of 8 kpc from 
the center is taken to be 5 H cm" 3 . The hydrogen 
density distribution is expressed by 

riH = cm cxp(— r/rn — z/zH)cos k (9 — rylog r/r a ) (6). 

Here, and r\ are the azimuthal angle and the pitch 
angle of the arms. 

4 Results of Simulations 

4.1 NGC 253 

The expansion velocity of the outer halo is estimated 
to be 200 km s _1 , which is required for shock heating 
of the gas up to 0.2 keV (2.3 x 10 6 K) (Pietsch ct al. 
2000). Taking this expansion velocity for the BHS at 5 
kpc height from the galactic plane in the ROSAT 0.75 
keV contour map, we find that the observation is well 
fitted by a shock front at ~ 3 x 10 7 yrs for an explosion 
energy of ~ 4 x 10 54 ergs as in Table 1. 

Among the models, however, a better fit is obtained 
by Models S2 and S6 than by the other models, which 
indicates that the decrease of halo gas density in the 
direction perpendicular to the disk is not as steep as 
in the other models. In Fig. 2 we show the shock 
front for Models S2 and S6 superposed on the 0.75 
keV ROSAT map of NGC 253 in contours overlaid on 
a B-band optical photograph from DSS. The ROSAT 
result shows that the hypershells tend to break and 
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become round at high altitudes, when the shock front 
reaches a few kpc from the galactic plane, mimicking 
a dumbbell shape. Hence, Model S6 appears to be 
preferable to S2, for which the BHS is more like an 
"hourglass" . Model S6 represents a case in which the 
gas density is uniform within the disk and halo up to 
a certain height (e.g. 2 kpc), where the halo faces a 
uniform intcrgalactic gas of lower density. Such density 
distribution will be possible if the halo gas is confined 
by the intergalactic gas pressure, but not stratified by 
a hydrostatic equilibrium under the disk gravity. 
— Fig. 2 — 

Although Model S6 appears better, overall reproduc- 
tion of the observations is obtained by both models. 
For simplicity, we hereafter use the Model S2 in which 
the shape of the shell can be well approximated by a 
part of an ellipsoidal sphere. This allows us a precise 
prediction for the X-ray emissivity and density distri- 
butions. Fig. 3 shows the calculated intensity distribu- 
tions for NGC 253 for two different opacity (k) values, 
which correspond to (a) a case of 0.75 keV band, and 
(b) 1.5 keV band. We also show the ROSAT X-ray im- 
ages in the right panels: right-top for 0.75 keV in gray 
scale, and the right-middle and right-bottom are for 
0.75 and 1.5 keV, respectively, overlaid on an optical 
photograph. As seen in the figures, the BHS simula- 
tion reproduces the following observed characteristics: 
(1) double-horn morphology of the diffuse emission fea- 
tures in the 0.75 keV band emerging from the disk, (2) 
strong absorption along the galactic plane, and (3) en- 
ergy dependence of the intensity distribution, namely, 
harder emission near the disk. 

- Fig. 3 - 

The Combination of three NGC 253 images in in- 
dependent energy bands centered at 0.25, 0.75 and 
1.5 keV enabled the construction of a true color im- 
age (Pietsch et al. 2000). For this purpose, the dif- 
ferent subbands were associated with the colors red, 
green and blue, for the 0.25, 0.75 and 1.5 keV band, 
respectively. Fig. 4 shows the thus obtained colored 
intensity distribution compared to the observed spec- 
tral intensity distribution taken from Pietsch et al. 
(2000), where red stands for soft and blue for hard. The 
present BHS model can nicely reproduce the observed 
spectrum distribution, such as the harder spectra in 
the disk and spiral arms due to stronger absorption of 
softer X-rays. 

- Fig. 4. - 

4.2 BHS at Various Inclination: Com- 
parison with Other Galaxies 

We have thus shown that the extended X-ray features 
in the halo of NGC 253 are well reproduced using the 



BHS model, which was originally used to model the 
North Polar Spur NPS and associated spurs in the 
Milky Way halo (Sofue 2000). It is, therefore, inter- 
esting to see, whether the BHS model also applies to 
extended X-ray features observed in other galaxies ex- 
hibiting starbursts and outflows. Here, we compare the 
observed morphology with the outcomes of our theo- 
retical description. We calculated 0.75 keV intensity 
distributions for various inclination angles at 0° (face- 
on), 30°, 45°, 60°, 80°, and 90° (edge-on). Fig. 5 shows 
the calculated results, and they are compared in Fig. 6 
with observed examples for M83 (i = 24°), NGC 1808 
(i = 58°), NGC 4258 (i = 67°), NGC 253 (i = 78°), 
NGC 3079 (edge-on), M82 (edge-on) and the Milky 
Way (edge-on). 

- Fig. 5 - 

- Fig. 6 - 

M83: When a galaxy with BHS is observed face- 
on, the X-ray images exhibit a bright ring feature as 
in Fig. 5. If the galaxy is slightly inclined such as 
i = 30°, the BHS ridges in both sides of the galactic 
plane are superposed to exhibit a box-shaped bright 
region. Such a box-shaped X-ray distribution appears 
to be observed indeed in the ROSAT PSPC image of 
M83, as shown in Fig. 6 (Strickland et al. 1997), where 
the arrow points toward the direction of the minor axis 
as the near side. Although it is difficult to distinguish 
the halo emission from disk contributions, spiral arms 
and/or the bar, some emission along the boxy edges 
could be due to overlapped ridges of BHS. 

NGC 1808: The extended 0.5-2 keV features bifur- 
cated in the halo of NGC 1808 as observed with ASCA 
(Iwaki et al. 1996) are also reproduced by our simula- 
tions, although the lopsidedness cannot be reproduced 
well. The arrow in Fig. 6 indicates the direction of the 
dusty jet along the minor axis in the nearer side. 

NGC 4^58: This is a galaxy known for its anoma- 
lous radio arms, along which X-ray extended emission 
is observed, as is shown in the 0.1-0.4 keV ROSAT map 
(Vogler and Pietsch 1999b). Although the resolution 
is still too crude, some spur features are mimicked by 
the BHS model. The origin of the anomalous arms is 
controversial, as described in section 2, but we com- 
ment that the anomalous radio arms (van Albada and 
van der Hulst 1982) are nicely fitted by two oppositely 
extending ridges of the BHS. In particular, their round 
shapes and curvature of the arms returning to the disk 
can be explained if they are part of BHS. However, the 
lobsidedness of the anomalous arms with respect to the 
rotation axis remains a question in the BHS model, al- 
though some lopsidedness is seen in the BHS of the 
Milky Way (Fig. 6 bottom, constructed from the en- 
ergy bands presented in Snowden et al. 1997). 

M82: The ROSAT HRI image (Strickland et al. 
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1997) shows a bipolar extended halo elongated in the 
direction perpendicular to the disk, and the emission 
appears to be centrally filled, indicating that the out- 
flow is well collimated, as well simulated by the "jet-like 
wind" model by Suchkov et al. (1994, 1996). We men- 
tion, however, that the northern end of the flow at a 
height of about 2.5 kpc (3') from the disk appears to 
be more opened into the intergalactic space, suggesting 
double-horn ridges or a partial shell. 

NGC 3079: The radio-lobe edge-on galaxy NGC 
3079 shows a box-shaped X-ray features, suggesting 
a double-horn morphology, as observed in the ROSAT 
PSPC image (Pietsch et al. 1998), and is well repro- 
duced by an edge-on BHS model. 

Milky Way: As mentioned in Sect. 2.6, the RASS 
reveals bipolar hypershells, which make up the North 
Polar Spur and its counter-spurs. The observations are 
well fitted by the BHS model. A detailed modeling for 
the case of the Milky Way is described in Sofue (2000). 

5 Discussion 

5.1 Probing Starburst and Halo Gas 
Distribution 

In the present shock approximation, the expansion ve- 
locity v at radius r is related to the mean ambient den- 
sity po and the explosion energy Eo as Eq ~ por 3 0v 2 . 
Here, r is measured from the observed shell radius, 
and the expansion velocity is also measured from the 
gaseous temperature inferred from X-ray spectrum. 
For given radii and velocities (temperature), which are 
both known parameters from observations, a scaling 
law can be established: In the case of constant E/p 
values, we obtain identical shock envelopes. 

For a given temperature (velocity), the emission 
measure is estimated from the intensity, and therefore, 
the accumulated total amount of gas in the shell is 
known, which is related to the mean gas density po- 
Hence, we can measure po, r and v from observations, 
and thus, we can finally estimate the total energy of 
the explosion. 

The morphology of the BHS is sensitive to the den- 
sity distribution in the disk and halo. The shell shape 
in the uppermost part manifests the halo-intergalactic 
density structure. If the intergalactic gas density is 
very low, the shell will have a larger on-average radius 
and a more conical form. In the case of higher den- 
sity, the shell becomes compact and round. Among 
the various ambient density distributions in Table 1, 
Models S2 and S6 appear to better fit the observations 
than the other models. Although both models allow 
a sufficient global reproduction of the observed mor- 
phology of NGC 253, Model S6 appears to be prefer- 



able to reproduce the dumbbell shape of the observed 
X-ray lobes. This suggests that the disk-to-halo gas 
density is rather uniform until it faces the intergalactic 
low-density region. Such uniform density distribution 
may be possible if the gas is confined by the gas pres- 
sure, not hydrostatic equilibrium in the disk gravity. 
In order to clarify this and fit the observations in more 
detail, higher resolution imaging in the soft X-rays is 
necessary. 

5.2 Origin of the Hypershells 

In our model, the energy injection at the galactic cen- 
ter is assumed to be impulsive, and the shock is strong 
enough to create a well-defined shell structure. The re- 
quired total energy given to the interstellar gas is of the 
order of ~ 10 55 ergs. The time scale of the explosive 
event responsible for the energy release is assumed to 
be significantly shorter than the expansion timescale of 
the hypcrshell, t ~ r/v ~ 10 7 years for r ~ several kpc 
and v ~ 300 km s^ 1 . So, the explosion (starburst) time 
scale must be shorter than a few million years. Such 
impulsiveness and robustness of the energy release can 
be explained if the nucleus of the galaxy has experi- 
enced a starburst ~ 10 7 years ago, lasting for a few 
million years or shorter, during which ~ 10 4 " 5 type II 
supernovae exploded. It is worth mentioning that sim- 
ilar timescales are predicted from starburst modeling 
based on infrared observations (e.g., Forster-Schreiber 
et al. 2000). 

Alternative models have also been discussed by sev- 
eral authors: Suchkov et al. (1994, 1996) developed a 
mass-loading outflow model, and simulated the X-ray 
jets in M82 using their hydrodynamical code. In their 
model, the blowout itself is as hot or as highly ener- 
getic as one would detect in the ROSAT band. This 
model may provide a better fit in the case of M82, 
where we observe a centrally-filled high energy flow 
well collimated in the direction perpendicular to the 
disk. A stellar-wind driven Galactic wind (Heckman 
et al. 1990) could also produce a shell structure, as ob- 
served with ROSAT for NGC 253 (Pietsch et al. 2000). 
If the duration of the starburst is short enough com- 
pared to the expansion time scale of the shell, namely, 
if the starburst is "impulsive" compared to the shell's 
life time, the wind model and the our shocked BHS 
model give essentially the same result. On the other 
hand, if the wind is steady and longer-lived than the 
shell's expansion time, e.g. a few tens of million years, 
the flow would become an open-cone shape, which does 
not apply to round shells observed in the case of NGC 
253. 
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6 Summary and conclusions 

We presented a theoretical model of the million degrees 
phase of the interstellar medium in the bulge, disk and 
halo of spiral galaxies. In this model, we assume a sud- 
den energy release in the nuclear region of the galaxy, 
as provided, e.g., by a nuclear starburst. The energy 
release is supposed to form a bipolar hyper-shell (BHS) 
of ionized ISM, which is heated to a million degrees and 
is thus visible in the soft X-ray band. Our simulations 
are based on the assumption of an adiabatic shock en- 
velope, and our projections of the three-dimensional 
data cubes for different galaxy inclinations take into 
account the absorption of soft X-rays by intervening 
cold gas in the disk of the galaxy. The model predicts 
different spectral energy distributions of the soft X-ray 
emission for different distances from the galactic plane. 

We applied our model to a mini-sample of nearby spi- 
ral galaxies, which are all classified as starburst galax- 
ies, and the model-predictions arc nicely supported by 
the observations. In the case of the very nearby galaxy 
NGC 253 (d = 2.6 Mpc), the X-ray images allow a high 
spatial resolution. Our model reproduces the ROSAT 
observations, namely the hourglass-shaped morphol- 
ogy, the absorption of disk and halo emission from the 
reverse (NW) side of the nearly edge-on galaxy due to 
cold gas in the disk, as well as the spectral distribution 
of the emission components. 

We extended our model to other spiral galaxies 
seen under different inclinations, edge-on (NGC 3079, 
M 82), under high inclination (NGC 4258, NGC 1808) 
or face-on (M 83). As for NGC 253, our simulations 
match the observed X-ray morphology and spectral en- 
ergy distribution. The X-ray corona of the Milky Way, 
observed in the ROSAT All Sky Survey, can also be 
matched with our model. This topic is briefly discussed 
in this publication and is subject of another, more de- 
tailed publication (Sofue 2000). 

From the comparison of our simulations with the ob- 
servations, we are able to deduce expansion velocities 
of the shells between 200 km s _1 and 300 km s _1 , which 
are necessary to heat the gas to (2 — 4) x 10 6 K. The age 
of the shells has to be of the order of some 10 7 years, 
and from this results we can derive an upper limit for 
the duration and intensity of the initial starburst which 
provided the energy release to form the BHS morphol- 
ogy. In general, this energy release is estimated to be 
of the order of 10 55 ergs for all galaxies, and typical 
timescales of the burst seem to be ~ several 10 6 years. 
Such timescales are also predicted from infrared obser- 
vations in the case of M 82, NGC 253 and NGC 1808 
(e.g., Forster-Schreiber et al. 2000). 

The new generation X-ray observatories, Chandra 
and XMM, both now in orbit, have larger energy bands 



than ROSAT had (0.2-10 keV versus 0.1-2.4 keV) and 
offer superior spectral resolution. They allow us to 
study spiral galaxies at sub-arcsec resolution ( Chandra) 
and to probe the ISM with an extremely high photon 
collecting power at a resolution of ~ 10" (XMM). The 
combination of Chandra and XMM results will provide 
spatially highly resolved images of the bright ISM fea- 
tures, as well as detailed X-ray spectra of the ISM in 
general. The BHS model, in which the shell morphol- 
ogy is sensitive to the ambient density distribution, will 
then allow us to probe the distributions of the million 
degree interstellar medium in the disk, halo and inter- 
galactic space as well as the disk halo interface at high 
spatial and spectral resolution. 

We are grateful to W. Pietsch for helpful comments 
and providing the NGC 3079 and NGC 253 data. We 
would like to thank Sue Madden for carefully proof- 
reading the manuscript. 
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Figure Captions 



Fig. 1: Calculated shock front in the galactic halo 
at every 2.5 x 10 7 yrs after an explosion and/or a 
starburst at the nucleus with a total energy of 0.4 x 
10 55 erg for various parameters as models SI to S6 
in Table 1 

Fig. 2: Calculated shock fronts for Model S2 and S6 
overlaid on the ROSAT 0.75 keV contours and op- 
tical photograph of NGC 253 from Pietsch et al. 
(2000) 

Fig. 3: Simulated intensity distributions for NGC 
253 at 0.75 and 1.5 keV (left panel) compared with 
ROSAT images at 0.75 and 1.5 keV taken from 
Pietsch et al. (2000) (right panel) 

Fig. 4: Addition of 0.25 (R), 0.75 (G) and 1.5 (B) keV 
intensities, compared with that of NGC 253 (Pietsch 
et al. 2000) 

Fig. 5: Simulated 0.75 keV views of a galaxy at various 
inclination angles 

Fig. 6: Comparison of the model with observed X-ray 
images of several galaxies with various inclinations: 
From top to bottom, M83, NGC 1808, NGC 253, 
NGC 4258, M82, NGC 3079, and the Milky Way. 
The arrows in M83 and NGC 1808 point the direc- 
tions of the minor axis in the near side. See the text 
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